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Summary 

The  temperatures  due  to  frictional  heating  within  a  solid 
lubricated  or  coated  journal  bearing  were  analyzed  by  using 
a  finite  element  method.  A  solid  model  of  the  shaft-bush 
tribocontact  was  generated  with  an  eight-node,  three- 
dimensional.  first-order  isoparametric  heat-transfer  element 
and  the  Patran  solid  modeler  software.  The  Patmar  (Patran- 
Marc)  translator  was  used  to  help  develop  the  Marc-based 
finite  element  program  for  the  system:  this  software  was  used 
on  the  Cray  X-MP  super  computer  to  perform  a  finite  element 
analysis  of  the  contact.  The  analysis  was  performed  for  various 
liner  materials,  for  thin.  hard,  wear-resistant  coated  bearings, 
and  for  different  geometries  and  thermal  cooling  boundary 
conditions.  The  analyses  indicated  that  thermal  conductivity 
of  the  liner  or  coating  material  is  the  most  vital  thermal 
parameter  that  controls  the  interf’ce  temperature.  In  addition 
to  design  variations,  the  proximity  of  the  cooling  source  to 
the  heat-flux-generating  interface  is  critically  important  to  the 
temperature  control  in  the  system. 

Introduction 

The  temperature  rise  due  to  frictional  heating  at  the  interface 
of  tribocontacts  has  attracted  the  attention  of  engineers  and 
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tribologists  in  recent  years.  Friction  and  w  ear  of  dry  contacts 
are  dependent  on  the  overall  temperature  increase  at  the 
interface.  Thermomechanical  effects,  such  as  the  initiation  of 
thermal  cracks  and  the  low-cycle  thermal-fatigue  failure  in  dry 
bearings  and  seals,  are  dependent  on  the  heat  dissipation 
mechanism  and  the  overall  temperature  rise.  Floquet.  Play, 
and  Godet  (ref.  1 )  used  a  two-dimensional  Fourier  transform 
method  developed  by  Ling  (ref.  2)  to  calculate  the  contact 
temperatures  in  dry  bearings  operating  with  a  plastic  liner. 
Later.  Floquet  and  Play  extended  this  technique  to  three- 
dimensional  problems  (ref.  3).  They  showed  that  small  changes 
in  the  design  may  bring  large  differences  in  contact 
temperatures,  and  that  to  reduce  the  maximum  temperature, 
the  rubbing  surfaces  must  be  as  close  as  possible  to  a  controlled 
temperature  heat  sink  and  must  not  be  shielded  by  any  thermal 
barrier.  Gecim  and  Winer  (refs.  4  and  5)  extended  the  Fourier 
transform  method  to  calculate  steady  temperatures  in  a  rotating 
cylinder  that  had  multiple  surface  heat  sources  and  a  variety 
of  geometries  and  thermal  boundary  conditions. 

In  recent  years,  high-temperature  applications  brought  about 
the  need  for  hard,  wear-resistant  coatings  (viz..  AUOj.  and 
ZrO;)  to  be  used  as  thermal  barriers  in  mechanical 
components.  When  these  components  are  subjected  to  heavily 
loaded  sliding  conditions,  "heat  checking"  often  occurs:  that 
is.  cracking  in  the  neighborhood  of  the  contact  zone  due  to 
the  combination  of  thermal  heating  and  mechanical  load.  This 
often  occurs  in  such  mechanical  devices  as  seals  and  brakes. 
Kennedy,  ct  al.  (refs.  6  to  8)  have  done  thermomechanical 
analyses  of  dry  sliding  contacts,  using  a  finite  element 
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technique.  In  particular,  a  recent  investigation  (ref.  8) 
addressed  the  problems  of  hard,  wear- resistant  coatings  on  a 
ductile  metallic  substrate,  us  applied  to  mecnanical  face  seals. 
A  thcrmo-elasto-plastic  analysis  of  stresses  and  deformations 
was  used. 

Ju  and  Chen  (ref.  9)  used  a  Fourier  transform  method  to 
investigate  the  mechanical  stresses,  temperature  field,  and 
thermal  stresses  due  to  a  moving  frictional  load.  Ju  and  Liu 
(refs.  10  and  11)  extended  this  technique  to  study  the 
phenomenon  of  thermomechanical  cracking  in  layered  media. 
Leroy.  Floquet.  and  Villechaise  (ref.  12)  also  used  a  fast 
Fourier  transform  algorithm  to  study  thermomechanical 
behavior  of  layered  media. 

Fourier  transform  techniques  generally  provide  useful 
analytical  expressions,  but  they  are  limited  in  applications 
dealing  w  ith  irregular  geometries  and  boundaries.  Insofar  as 
we  know,  dry  bearing  temperatures  have  not  been  investigated 
with  the  finite  element  technique  that  accommodates  irregular 
boundaries.  Herein,  results  of  a  finite  element  analysis  of  dry 
contact  temperatures  due  to  frictional  heating  at  the  interface 
of  a  shaft-bush-housing  system  are  presented.  A  solid  model 
of  the  system  was  generated  with  the  Patran  solid  modeler 
software,  and  the  finite  element  analysis  was  done  with  the 
Marc  finite  element  software.  The  temperature  distribution  in 
the  system  was  determined  with  various  liner  materials,  coating 
materials,  coating  thicknesses,  and  bearing  geometries. 

This  work  also  shows  that  the  method  can  serve  as  a 
computer-aided  design  technique  to  model  mechanical 
components  (such  as  the  cylinder-piston  of  an  advanced  low 
heat  rejection  (LHR)  diesel  engine  or  Stirling  engine,  and 
aircraft  control  and  undercarriage  bearings)  for  temperature 
rise  and  thermal  distortions  under  specific  operating 
conditions. 


Symbols 

IB]  temperature  gradient  interpolation  matrix 

[C|  capacitance  matrix 

c  specific  heat.  J/kg  °C 

/  friction  coefficient  at  the  contact 

li  film  heat  transfer  coefficient.  W/nr  °C 
|K|  heat  transfer  or  coefficient  matrix 

k  thermal  conductivity.  W/m  °C 

M  number  of  finite  elements 

[N|  temperature  interpolation  matrix 

n  unit  vector 

p  contact  pressure.  Pa 

Q  internal  heat  generation 

q  heat  tlux  distribution.  W/nr  or  N/cm-scc 


Rb,  inner  radius  of  the  bearing,  m 
>-  number  of  nodes  in  an  element 

S  surface  area,  m- 

T  temperature.  °C 

/  time,  sec 

V  rubbing  velocity,  m/sec 

x.y.z  Cartesian  coordinates 
T  surface  of  integration,  nr 

k  thermal  diffusivity.  nr/sec 

p  density,  kg/m 3 

Q  volume  of  integration,  nr' 

uy  angular  speed  of  the  shaft,  rad/sec 

Subscripts: 

c  conductance 

c  element,  in  a  global  sense 

li  surface  convection 

m  maximum 

q  surface  heating 

.v  surface 

T  nodal  temperatures 

Superscript: 

c  element,  in  a  local  sense 

Theory 

The  theory  discussed  herein  is  taken  from  reference  13. 

Finite  Element  Formulation 

The  steady-state  or  transient  heat  transfer  in  a  three- 
dimensional  anisotropic  solid  moving  with  velocity  V  is 
governed  by  the  Fourier  heat  conduction  or  energy  equation, 
which  is  written  as 


v  -k  vT  -  pc 


1 y-vT  +  (?=() 


(1) 


where  Q(x,y,z,t)  is  the  internal  heat  generation  rate  per  unit 
volume.  The  solution  for  the  temperature  distribution  must 
satisfy  the  boundary  conditions  and  initial  conditions  on  all 
portions  of  the  surface.  The  inertia  conditions,  which  specify 
the  temperature  at  time  zero,  are  required  for  the  transient 
analysis.  The  heat  conduction  boundary  conditions  are  as  follows: 


where 


T=  7j 

qn  =  —  kvT-ii 
h  ( Ts  -  7",.)  n  =  k  v7-n 

where 

7j  prescribed  temperature  in  the  specified  region 

S,  of  the  surface 

—  K^T-n  prescribed  heat  flux  normal  to  the  surface  S2 
k  v  7'-n  convection  to  exchange  temperature  7j,  from  the 
surface  with  unknown  temperature  Ts 

The  solution  domain  in  the  solid  is  discretized  into  M  finite 
elements  possessing  r  nodes.  Temperature  and  temperature 
gradients  within  the  element  are  expressed  in  terms  of  nodal 
temperatures  as 


and 


Te(x.v,z,t)  =  lN(.v,v,z,r)]  T(t ) 


dT 

—  ix,y.z.t) 
dt 


dT 

—  (.r.y.c,/) 
dt 


=  [B(.v,y,;)]  T(t) 


(3) 


dT 

—  (r.v.c.r) 

dt 


where  [N]  is  the  temperature  interpolation  matrix,  and  [B]  is 
the  temperature  gradient  interpolation  matrix  such  that. 
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Finite  element  equations  are  derived  from  the  energy  equation 
(eq.  1 )  by  using  the  method  of  weighted  residuals  (Galerkin's 
method): 

[C]  j^j  +[|Kfl  +  [K*||  17] 

=  Ry  4-  R(/  +  R>,  +  R7-  (5) 
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The  coefficient  matrix  [C]  of  the  time  derivative  of  the  nodal 
temperatures  is  the  element  capacitance  matrix.  The  coefficient 
matrices  [K,  ]  and  [K,,]  are  element  conductance  matrices  and 
relate  to  conductance  and  convection. 

The  vectors  Rr .  R0 .  R(/ .  Rh  are  heat-load  vectors 
arising  from  specified  nodal  temperatures,  internal  heat 
generation,  specified  surface  heating,  and  surface  convection, 
respectively.  For  steady-state  analysis,  equation  (3)  reduces  to 


[[Kc]  +  IK*]]  T  =  Rq  +  P,  +  R„  (7) 


Problem  Statement  and  Solution  Methodology 


In  the  present  work  a  shaft-bush/liner-housing  type  of 
tribosystem  was  analyzed.  The  system,  shown  in  figure  1. 
consists  of  a  hollow  shaft  (1)  rotating  in  a  bush-liner  (2)  fitted 
within  a  stationary  housing  (3).  Alternatively,  the  shaft  rotates 
in  a  housing  with  a  thin.  hard,  wear-resistant  ceramic  coating 
at  its  internal  surface.  There  is  frictional  rubbing  at  the 
interface  of  the  liner  or  coating  and  the  shaft.  Frictional  heating 
at  the  interface  is  conducted  to  the  shaft  and  the  bearing  and 
housing,  thereby  causing  a  temperature  rise  in  both.  The  main 
objective  of  the  present  study  was  to  develop  a  finite  element 
model  of  the  previously  described  system  by  using  the  Patran 
solid  modeler  software,  and  thereafter  to  obtain  a  Marc  finite 
element  program  by  using  a  Patran-Marc  (Patmar)  translator. 
Following  this,  we  analyzed  the  problem  for  temperature 
distribution  within  the  solid,  using  the  Marc  finite  element 
software,  subject  to  the  proper  boundary  conditions  and  the 
appropriate  material  thermal  properties.  The  results  were 
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postprocesscil  with  the  Patran  postprocessor.  The  boundary 
conditions  tor  the  system  analyzed  were 

( 1)  A  heat-transfer  coefficient  h  equal  to  58  W/nr  °C  and 
a  reference  temperature  of  20  °C.  caused  by  natural 
convection 

(2)  An  imposed  temperature  of  35  °C  at  the  outer  ends  of 
the  shaft,  because  of  direct  cooling 

(3)  An  imposed  temperature  of  35  °C  at  the  periphery  of 
the  housing,  because  of  direct  cooling 

(4)  An  imposed  temperature  of  35  °C  for  Case  I.  caused 
by  direct  cooling  at  the  inner  radius  of  the  shaft 

Shaft  speeds  of  2.6.  6.5.  and  13.0  rad/sec  were  used  to 
calculate  the  heat  flux  generated  at  the  interface.  A  nodal  heat 
flux  q  equal  to  the  product  of  the  contact  pressure  p.  the  friction 
coefficient  /.  and  the  rubbing  velocity  V  w'as  input  to  the 
interface  between  the  shaft  and  bush  (or  coated)  housing.  An 
harmonic  contact  pressure  distribution  extending  from  ±90° 
was  assumed  (i.e..  the  journal  contacts  of  half  the  bearing). 

Effects  of  various  liner  materials,  coating  materials,  coating 
thicknesses,  and  geometries  on  the  temperature  distribution 
and  the  maximum  temperature  at  the  interface  were 
investigated. 

Results  and  Discussion 

The  maximum  interface  temperatures  and  the  temperature 
distribution  in  the  shaft-bearing  system  are  presented  herein. 
This  system  was  composed  of  either  a  bearing  with  liner 
material  in  a  housing,  or  a  bearing  with  a  thin.  hard,  wear- 
resistant  coating  on  its  inner  surface.  The  temperature  rise  was 
due  to  frictional  heating  at  the  interface.  The  hollow  shaft  had 
different  dimensions  (see  table  I)  for  the  three  cases  studied: 

( 1 )  bearings  with  soft  liner  materials:  (2)  bearings  with  thin, 
hard  coatings  on  the  internal  surface:  and  (3)  variations  in 
geometry  (table  II)  and  the  cooling  condition.  Nog  u  heat  flux 
generated  in  the  contact  was  estimated  to  be  equal  to  the 
product  of  the  contact  pressure,  rubbing  velocity,  and  friction 
coefficient.  The  coefficients  of  friction  for  the  liner  and  coating 
materials  ranged  from  0. 15  to  0.3.  However,  an  intermediate 
value  of  0.2  was  assumed.  Our  objective  was  to  evaluate  the 
temperature  rise  in  the  contact  for  the  same  amount  of  heat 
flux  input.  A  harmonic  pressure  distribution  extending  from 
±90°  was  assumed. 

Thus,  heat  flux  distribution 

q  =  q,„  cos  0 

where 

</»,  =  .fVPm 
and 

/  friction  coefficient  at  the  contact,  assumed  to  be  0.2 
V  rubbing  velocity:  Rhl u.\ 
p,„  maximum  contact  pressure 


Case  I— Bearings  with  Soft  Liner  Materials 

Thermal  properties  of  various  liner  materials  are  show  n  in 
table  111.  The  shaft  and  the  housing  were  made  of  steel.  For 
a  6.0-mm  thick  liner,  the  maximum  interface  temperatures 
obtained  for  various  liner  materials  are  shown  in  table  IV.  The 
lowest  of  the  maximum  temperatures  was  108.5  °C  for  the 
graphitized  metal  liner,  which  had  the  highest  thermal 
conductivity  and  diffusix ity  values.  As  the  material  s 
conductivity  and  diffusivity  decreased,  the  temperature  at  the 
interface  increased.  In  the  absence  of  conv  ection,  the  thermal 
conductivity  of  the  liner  was  the  material's  single  most 
important  thermal  property .  influencing  both  the  temperature 
rise  at  the  interface  and  the  temperature  distribution  as  well 
(fig.  2).  Table  V  shows  the  variation  of  maximum  interface 
temperature  with  five  different  heat  distribution  functions  for 
a  resin  liner  material  (ref.  1 ).  w  hile  keeping  all  other  conditions 
constant.  We  observed  that  the  maximum  interface  temperature 
and  the  temperature  distribution  can  differ  significantly  w  ith 
the  nature  of  the  heat  distribution  in  the  contact  region  (fig.  3). 
The  nature  of  the  heat  distribution  depends  on  the  contact 
pressure  in  the  contact  region,  which  in  turn  depends  on  the 
geometry  and  the  type  of  loading.  The  nature  of  the  loading 
may  vary  with  the  operating  conditions. 

The  effects  of  the  magnitude  of  the  load  and  the  thickness 
of  the  liner  material  are  show  n  in  table  VI.  Not  surprisingly  . 
a  thicker  liner  material  can  result  in  a  lower  interface 
temperature,  provided  the  conductivity  of  the  liner  material 
is  higher  than  the  conductivity  of  the  bush  material.  Thus,  for 
higher  loads  or  pV  factors,  a  thicker  graphitized  metal  liner 
would  help  reduce  the  temperatures. 

Case  II— Bearings  with  a  Thin.  Hard  Coating  on  the 
Internal  Surface 

Thermal  properties  of  thin.  hard,  coating  materials  are  giv  en 
in  table  VII.  For  a  1.0-mm-thick  coating,  the  variations  in 
maximum  interface  temperatures  are  show  n  in  table  VIII.  As 
is  evident  from  table  VII  for  the  thermal  properties  and  table 
VIII  for  the  maximum  temperatures,  the  thermal  conductiv  ity 
of  the  coating  material  determines  the  magnitude  of  the 
maximum  interface  temperature.  When  a  coating  material  with 
low  thermal  conductivity  was  used  as  a  thermal  barrier  to 
shield  the  bearing,  high  temperatures  occurred  at  the  shaft¬ 
bearing  interface,  thereby  dissipating  more  heat  to  the  shaft 
and  causing  higher  shaft  temperatures. 

In  order  to  study  the  effect  of  coating  thickness,  calculations 
were  made  for  titanium  nitride  coatings  ranging  from 
0.125  to  1.0  mm  thick.  The  resulting  maximum  interface 
temperatures  arc  shown  in  table  IX.  Note  that  maximum 
interface  temperatures  decreased  with  decreasing  coating 
thickness.  This  means  that  thicker  coatings  with  lower 
conductivity  than  the  substrate  (such  as  titanium  nitride  on 
steel)  are  better  thermal  barriers  than  very  thin  coatings. 
Temperature  distributions  in  the  system,  which  result  from 
usin>’  three  different  thicknesses  of  titanium  nitride  coating. 
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are  shown  in  figures  4  to  6.  Similar  calculations  for  various 
thicknesses  of  tungsten  carbide  coatings  showed  almost  no 
change  in  maximum  interface  temperature,  since  tungsten 
carbide  has  a  higher  thermal  conductivity  than  the  steel 
substrate. 

Case  III— Variations  in  Geometry  and  the  Cooling 
Condition 

Tables  X  and  XI  show  the  variation  of  maximum  interface 
temperatures  for  a  number  of  bearing  outer  radii  and  bearing 
lengths.  Other  dimensions  remained  the  same.  Higher  interface 
temperatures  occurred  with  larger  bearing  diameters,  since 
the  cooling  source  was  located  at  the  outer  bearing  diameter. 
An  increase  in  bearing  length  did  not  alter  the  maximum 
temperature.  For  a  hollow  shaft  that  was  cooled  internally, 
the  interface  temperature  was  drastically  reduced— from  203 
to  109  °C.  The  thermal  properties  of  the  materials  are 
considered  to  be  temperature-independent  here.  Figures  7  to 
10  show  the  changes  in  the  temperature  distribution  in  the 
system  because  of  geometry  variation  and  inner  surface  cooling 
of  the  hollow  shaft. 

An  increase  in  rubbing  velocity  or  pV  factor  can  cause  a 
tremendous  rise  in  the  maximum  interface  temperatures.  Table 
XII  shows  that  temperatures  as  high  as  928.0  °C  can  be 
obtained  at  pV  =  5362.5  N/cm-sec  with  a  1.0-mm  thick 
titanium  nitride  coating.  Figure  11  illustrates  the  thermal 
distribution  for  this  case. 


Summary  of  Results 

A  Patran  solid  modeling  package  was  used  to  generate  a 
finite  element  model  for  the  thermal  analysis  of  a  shaft-bush 
tribosystem.  The  Marc  finite  element  software  was  used  to 
determine  the  temperature  distribution  in  the  tribosystem. 
These  software  permit  finite  element  modeling  of  bearings  and 
seals  so  that  temperature  distributions  can  be  calculated  in 
specific  applications.  Mechanical  components  (viz.,  the 
cylinder-piston  of  an  advanced  LHR  diesel  engine,  the  Stirling 
cycle  engine  for  space  powei  applications,  aircraft  control  and 
undercarriage  bearings,  and  so  on)  can  be  modeled  for 
temperature-rise  and  thermal  distortions  under  specific 
operating  conditions.  The  results  of  our  analysis  can  be 
summarized  as  follows: 

1.  If  a  constant  friction  coefficient  is  assumed,  thermal 
conductivity  of  the  liner  or  coating  material  is  the  most 
important  thermal  property  for  controlling  temperatures  in 
tribocontacts. 

2.  Design  changes  in  bearing  diameter  and  length  can  result 
in  significant  differences  in  the  interface  temperature. 

3.  Interface  temperatures  can  be  substantially  reduced  by 
cooling  the  inner  surface  of  the  hollow  shaft. 

4.  Designs  that  can  provide  a  cooling  medium  in  close 
proximity  to  the  heat-flux-generating  interface  would  be  most 
effective  in  reducing  the  temperature  rise  in  the  contact. 
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TABLE  I, -SHAFT  SPECIFICATION  AND  BEARING 
DIMENSIONS  FOR  THREE  CASES 


Case 

Shaft 

Bearing 

Outer 

Inner 

Length, 

Outer 

Inner 

Length. 

radius. 

radius. 

cm 

radius. 

radius. 

cm 

cm 

cm 

cm 

cm 

aI  and  JII 

2.5 

1.5 

8.4 

6.0 

2.5 

1.2 

hIII 

2.5 

1.5 

16.0 

5.0 

2.5 

1.6 

_ 

is  sot  oil  shall  to  ri$:hl  hs  I  2  cm 
^Bearing  is  ccnlraflv  liK'atcd  on  shall. 


TABLE  II.— GEOMETRY 
VARIATIONS  FOR 
CASE  III 


TABLE  IV. -MAXIMUM 
INTERFACE  TEMPERATURES 
FOR  VARIOUS  LINER 
MATERIALS 


[Heal  flux  distribution  q  =  q„ cosd>: 
pV  =  1072.5  N/cm-spc;  liner 
thickness  =  6.0  mm.) 


Liner  material 

Maximum 

temperature. 

°C 

Graphitized  metal 

108.5 

PTFE-bronze 

118.6 

Carbon  graphite 

147.0 

PTFE-thermoset 

152.8 

PTFE-ceramic 

152.7 

MocS-nylon 

152.8 

Resin 

_ 

152.8 

Bearing 

Bearing  outer 

length. 

radius. 

cm 

cm 

1.6 

3.2 

1.6 

5.0 

16 

8.0 

3.2 

5.0 

TABLE  V. -MAXIMUM  INTERFACE 
TEMPERATURES  FOR  VARIOUS  HEAT 
FLUX  DISTRIBUTION  FUNCTIONS 

[Resin  liner  thickness  =  6.0  mm:  pV  =  1072.5 
N/em-sec.) 


TABLE  III. -THERMAL  PROPERTIES  OF  LINER  MATERIALS 


Heat  flux  distribution. 

<7- 

N/cm-sec 

Maximum 

temperature. 

”C 

Maximum 

temperature 

location. 

deg 

</mcos  <S> 

152.8 

0 

<7„,-(4/7r)cos:d> 

179.7 

0 

c/,„  *(2/70(1  +sin</>) 

169.7 

±90 

qm’(2/ir) 

114.0 

0 

qm'  (2/tt)(  1  -sin<(>) 

99.6 

0 

|From  ref.  14.) 


Liner  Material 

Specific 

heat. 

c. 

J/kg  °C 

Thermal 

conductivity, 

k. 

W/m  °C 

Density. 

P- 

kg/m-1 

Diffusivity. 

KlpC. 

m:/sec 

Graphitized  metal 

377.0 

127.4 

7039  0 

4.8  x  10  ^ 

PTFE-bronz.c 

394.0 

42.5 

6501.6 

1.66  x  10  5 

Carbon  graphite 

754.0 

3.35 

2101.7 

2.1  x  10" 

PTFE-thermoset 

942.8 

.263 

2293.2 

1.1  x  IO-fi 

PTFE-ceramic 

838.0 

.337 

2376.0 

1.69  x  10-7 

Mo^S-nylon 

188.6 

.25 

1198.0 

1.22  x  10-7 

Resin 

1214  0 

.233 

1200.0 

1.6  x  10 -7 

TABLE  VI. -MAXIMUM  INTERFACE  TEMPERATURES 
FOR  pV  FACTORS  AND  LINER  THICKNESSES 

|Heat  flux  distribution  q  =  qm  cos  <2>. ] 

(b)  Graphitized  metal  liner; 

(a)  Resin  liner.  6.0  mm  thick  pV  =  1072.5  N/cm-sec 


Liner 

Temperature. 

thickness. 

°C 

mm 

;  I 

1 

18.0 

ptrlV  factor. 

Temperature. 

N/cm-sec 

°C 

1072.50 

152.8 

536.25 

93.9 

107.25 

46  8 

6 


TABLE  VII  -THERMAL  PROPERTIES  OF  COATING 
MATERIALS  AND  SUBSTRATE 


Material 

Specific 

heat. 

c. 

J/kg  °C 

Thermal 

conductivity. 

k. 

W/m  °C 

Density. 

P- 

kg/m1 

Diffusivity 

KlpC. 

m-/scc 

Aluminum 

869.8 

267.7 

2  699.0 

1.I36X  10  1 

Zirconium 

17  769.0 

131.9 

6  513.0 

1.142x10  h 

Silicon  carbide 

711.0 

104.4 

2  997.0 

4.9x  10  s 

Tungsten  carbide 

100.0 

100.0 

15  800.0 

5.0x10  6 

Aluminum  oxide 

1  039.0 

20.73 

3  997.0 

4.99 X  I0-o 

Stellite 

421.9 

9.7 

8  300.0 

2.77x10  -« 

Titanium  nitride 

590.0 

4.0 

5  200.0 

1.3xl0-6 

Zirconia“ 

400.0 

1.8 

5  750.0 

7.8x10  s 

Steel 

460.0 

46.0 

7  800.0 

l.28x  I0-* 

uParnall>  stabilized 


TABLE  VIII. -MAXIMUM 
INTERFACE  TEMPERATURES 
FOR  COATING  MATERIALS 


[Coating  thickness  =  1 .0  mm;  heal 
flux  distribution  q  =  qmcos<t>: 
pV  =  1072.5  N/cm-sec.) 


Coating  material 

Maximum 

temperature. 

”C 

Aluminum 

176.8 

Zirconium 

181.1 

Silicon  carbide 

182.2 

Tungsten  carbide 

182.4 

Aluminum  oxide 

191.3 

Stellite 

199.5 

Titanium  nitride 

213.0 

Zirconiaa 

234.7 

TABLE  IX. -MAXIMUM 
INTERFACE 
TEMPERATURES  FOR 
VARIOUS  TITANIUM 
NITRIDE  COATING 
THICKNESSES 


[Heat  flux  distribution 
q  -  q,„cos<t>;  pV 
=  1072.5  N/cm-sec. j 


Coating 

thickness, 

mm 

Maximum 

temperature. 

°C 

1.0 

213.0 

.75 

208.0 

.5 

201.0 

.25 

194.0 

.125 

190.0 

aPar* liillv  Mahdi/cd 


TABLE  X. -MAXIMUM  INTERFACE 
TEMPERATURES  FOR  VARIOUS 
BEARING  DIMENSIONS 

|Coating.  1-mm-thick  TiN:  heat 
(lux  distribution  q  =  q,„ cosip; 
pV  =  1072.5  N/cm-sec.  1 


Bearing  and  housing 

Maximum 

outer  radius. 

temperature. 

cm 

°c 

3.8 

182.0 

5.0 

202.1 

8.0 

228.1 

TABLE  XL— MAXIMUM 
INTERFACE  TEMPERATURES 
FOR  VARIOUS  BEARING 
LENGTHS  AND  COOLING 
CONDITIONS 


(Coating,  l-mm-thick  TiN; 
pV  =  1072.5  N/cm-sec;  bearing 
outer  radius.  5  cm.) 


Bearing  length. 

Maximum 

cm 

temperature. 

°c 

1.6 

202.1 

3.2 

202.5 

A3. 2 

109.0 

■“With  shaft  inner  surface  cooling. 


TABLE  XII. -MAXIMUM 
INTERFACE 
TEMPERATURES  FOR 
VARIOUS  pV  FACTORS 


[Coating.  1-mm-thick  TiN; 
heat  flux  distribution 
q  =  q„,  cos<f>.| 


pV  factor. 

Maximum 

N/cm-sec 

temperature. 

°C 

1072.5 

213.0 

2681.25 

481.7 

5362.5 

928.0 

TompoMtuM? 

C 


(teure  ’  I  rihosNsicin  tcinperaliire  ci'nlours  lor  case  I  Bearing:  has  h  0  mm  ihich  soli  resin  linei 


1'iiMitc  4  I'rthoswcni  lotnpor.ilurc  contours  lor  c.im:  II  Bc.irinj:  h.i>  I  O  mm  thick  li\  outing 


Temperature 

"C 


•Mi  U'M)|hi .it mi i'  ti*i  v.im*  111  |-ii*.tt in**  ti.is  I  n  mm  ihuk  li\  v-'.ifiru-  Iv.itm/  i.i.t 
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